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Abstract. We have surveyed an optical/IR selected sample of nearby E/SO galaxies with and without nuclear dust structures 
with the VLA at 3.6 cm to a sensitivity of 100 (jjy, We can construct a Radio Luminosity Function (RLF) of these galaxies 
to ~ 10 19 W Hz -1 and find that ~ 50% of these galaxies have AGNs at this level. The space density of these AGNs equals 
that of starburst galaxies at this luminosity. Several dust-free galaxies have low luminosity radio cores, and their RLF is not 
significantly less than that of the dusty galaxies. 
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1. Introduction 

Research conducted during the last decade gave a new look to 
nearby elliptical galaxies previously considered as old, uniform 
systems with little gas or dust. Images from the Hubble Space 
Telescope (HST) have shown that many early-type galaxies 
have a large amount of dust (10 3 — 10 7 M ), either in the 
form of a nuclear disk or in a more diverse shape of filaments. 
Among different studies there is a large variation in the detec- 
tion rates which may be due to the different methods, resolu- 
tions, and sensitivities of the observations (Sadler & Gerhard 
1985 40%; Goudfrooij et al. 1994 41%; van Dokkum & Franx 
1995 48%; Ferrari et al. 1999 75%; Tomita et al. 2000 56%; 
Rest et al. 2001 43%; Tran et al. 2001 (IRAS bright sample) 
78%), but the general conclusion is that dust is common in 
nearby ellipticals. 

Establishing the presence of dust in nearby early type 
galaxies is only the first step towards determining the role of 
dust in these systems. It is already a well-known fact that radio- 
loud ellipticals often have large amounts of dust but there are 
some open questions, specially for the radio weak sources. 
Verdoes Kleijn et al. (1999) found that the incidence of dust in 
radio-loud early type galaxies is 89% while Tran et al. (2001) 
has a value of 43% for the occurrence of dust in their snapshot 
sample of relatively radio-quiet nearby early-type galaxies (for 
a description of the sample see Sec. 2). In the same sample, 
66% of dusty galaxies have NRAO VLA Sky Survey (NVSS) 
1.4 GHz flux detections (Condon et al. 1998), while only 8% 
of galaxies without dust are listed as radio sources. 
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These results raise a question: how important is the pres- 
ence of dust for radio emission in the nuclei of ellipticals? 
Plausibly, dust indicates the presence of gas, and gas is nec- 
essary to fuel the activity of a central massive black hole (BH). 
However this line of reasoning is highly incomplete. Gas may 
be present without dust. Dust may be present but not visually 
detectable (Goudfrooij & de Jong 1995). Dust and gas that have 
fed a BH in the past may not be observable at the time when the 
nuclear activity is observed. These arguments justify a careful 
study of the relation between dust and nuclear radio emission to 
determine the relevance of radio luminosity, dust morphology 
and other effects. 

There are two approaches to the study of extragalactic ra- 
dio sources. The first one is based on catalogs of discrete ra- 
dio sources followed by an analysis of the optical counter- 
parts. The second involves searching for radio emission from 
optically chosen objects. The first approach (e.g. de Koff et 
al. 2000), is relatively efficient in finding radio galaxies, but 
emphasizes powerful radio sources and may not provide a good 
counter sample of radio quiet galaxies. The second approach 
conversely emphasizes weak radio sources (e.g. Sadler, Jenkins 
& Kotanyi 1989; Wrobel 1991; Wrobel & Heeschen 1991; 
Sadler et al. 2002). 

Both types of radio surveys are important. Here we have 
chosen the second method primarily so that the optical selec- 
tion of the sample, including Hubble type and especially dust 
content is not biased by a priori selection for radio emission or 
other "interesting" properties of the galaxies. The survey ob- 
jects are selected on their optical/IR properties only and then 
observed with the VLA with the purpose of establishing the 
presence of nuclear AGNs. We compare our dusty and non- 
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dusty parts of the sample to investigate the importance of dust 
(as a fuel reservoir) for the existence of nuclear activity. 

In Section 2 we present the sample and discuss the observa- 
tions and the data reduction. In Section 3 we present the results 
of our study. They are followed with a discussion in Section 4. 
Section 5 brings a discussion on correlation of dust with radio 
emission. The conclusions are given in Section 6. 

2. Observations 

2. 1. Sample 

Our sample is compiled from two different samples described 
by Rest et al. (2001) and Tran et al. (2001). The first sample was 
created by selecting E/SO galaxies on their optical properties 
only from the Lyon/Meudon Extragalactic Database (LEDA). 
A randomly selected subset of 68 galaxies from this sample 
were observed with HST using WFPC2 in snapshot mode and 
thus this sample is referred to as the "snapshot" sample. An ad- 
ditional sample of galaxies was assembled from archival HST 
images of nearby E/SO galaxies selected for their 100 /jm IRAS 
emission as these were likely to contain large amount of dust 
(Tran et al. 2001). This sample is referred to as the "IRAS sam- 
ple". From these two samples, we have selected 36 objects (18 
from each) according to their optical/IR properties, with no 
regard to radio properties. The reason for this selection was 
to avoid biasing in picking a priori "interesting" objects and 
objects with strong radio fluxes. Galaxies in our sample are 
nearby ellipticals and lenticulars (E/SO), cz < 3200 kms -1 , 
at galactic latitude exceeding 20° to minimize Galactic ex- 
tinction, and with absolute V-band magnitude less then -17. 
Because of their optical/IR selection they tend to have low ra- 
dio powers. 

The global properties of galaxies in our sample are listed in 
Table 1 of Rest et al. (2001) and in Table 7 of Tran et al. (2001). 
In the list of 36 galaxies, 18 of them were chosen because they 
have dust in the form of disks or filaments. The other 18 non- 
dusty galaxies were selected to match dusty galaxies in optical 
properties, redshift, magnitude, and IRAS flux. However, af- 
ter the initial selection, more detailed studies (Rest et al. 2001, 
Tran et al. 2001) showed that 6 of the "non-dusty" galaxies 
showed faint dust structures and have here been included in the 
"dust" class. We used H=80 kms -1 Mpc -1 to be consistent 
with the papers defining the samples. 

2.2. Data Acquisition and Reduction 

The observations were undertaken with the VLA in C con- 
figuration at 3.6 cm wavelength. All sources were observed 
at two frequencies in the 8 GHz X-band (8.4351 and 8.4851 
GHz) with a bandwidth of 50 MHz for each frequency. We ob- 
served 68 sources in total, 36 galaxies and 32 calibrators. Each 
galaxy was observed for 15 minutes while calibrators were ob- 
served for 130 seconds. Most of the calibrators had position 
code A (positional accuracy < 07002), but four calibrators had 
B (0 / .'002 - O'.'Ol) and three had C (0701 - 0'.'15) as is indicated 
on the calibrator web page of the VLA. The radio positions of 
the detected sources are limited by this positional accuracy of 



the calibrators, as well as by the accuracy of the Gaussian fit 
to the source brightness distribution, which is dependent on the 
signal-to-noise ratios. Taking this in account the overall accu- 
racy is about 50 mas for mJy sources and about 100 mas for 
100 fiiy sources. The observations were taken on March 13, 
2000. 

We used the Astronomical Image Processing System 
(AIPS) to reduce the data using the standard procedures from 
the AIPS cookbook. After initial calibration, the data were im- 
aged using the task IMGR. The data were self-calibrated in 
phases to improve the image dynamic range, using a model 
derived from the same data. In some cases amplitude self- 
calibration was performed on the data to improve the final im- 
ages. For our astrometric purpose, the positions of the sources 
were extracted before self-calibration so that phase information 
was preserved. All the images were examined using the tasks 
JMFIT and IMSTAT. 

3. Results 

Twenty galaxies in our sample of 36 were detected as radio 
sources. Three detected sources (associated with NGC 2986, 
NGC 3610, NGC 4125) can not be matched with the cen- 
tral regions of the galaxies and there are no visible counter- 
parts on the available HST pictures, hence they are most likely 
background sources. The radio sources lay far from the nu- 
clei (about 2f67for NGC 2986, 3.'84for NGC 3160, and l'.U 
for NGC 4125). Although the fluxes and positions of these 
sources are listed in Table 1 (with asterisks) we treat them as 
non-detections of central AGNs in the surveyed galaxies. This 
leaves 17 AGN detections in 36 galaxies (47 % detection rate). 
The smallest signal to noise ratio (SNR) is about lOcx with a 
survey average rms a = 2.8 x 10~ 5 Jy/beam. For non-detected 
sources we calculated 4cr upper limits on detection, thus, the 
detection limit of our survey is about 0.1 mJy. Radio properties 
of the sample are given in Table 1 . By comparison the detec- 
tion limit of the NVSS (Condon et al. 1998) used by Tran et 
al. (2001) to discuss radio properties of our sample is <~ 3 mJy, 
a factor of 30 higher. 

Most of the detections are point-like, unresolved structures. 
NGC 5322 is the only galaxy with noticeable jet-like struc- 
ture. Typical detected sources are on the level of a few mJy; 
the weakest detections were <~ 200^iJy. Of the 36 galaxies in 
the sample, 24 galaxies show disk or filamentary dust structure 
and 13 (54%) of them are detected as radio sources. Twelve 
show no dust of which four (33%) are detected. 

4. Discussion 

4. 1. Nature of Detected Radio Sources 

Most of the detections are unresolved radio sources easily as- 
sociated with the central 175 on the HST image. At 25 Mpc, 
the mean distance of the galaxies in the sample, 1" is about 
120 pc. Thus the emission is clearly (near) nuclear, but not 
necessarily of AGN origin. Since the sources are weak (radio 
power ranges from 10 19 W Hz -1 to 10 21 W Hz -1 with a few 
higher exceptions) there is a possibility that they arise from a 
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Table 1. Radio properties of the galaxies. 



name 


dust 


D 


peak flux 


RA 


DEC 


L 


S 


NVSS flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ngcl400 


2 


25.4 P 


2.092 ± 0.02 


03 39 30.815 


-18 41 17.42 


1.61 ± 0.02 


1.80 


2.5 ± 0.5 


ngcl439 


4 


20.9 T 


<0.1 












ngc2549 





15.7 B 


<0.1 












ngc2592 


4 


25.5 H 


0.41 ± 0.02 


08 27 08.040 


25 58 13.00 


0.32 ± 0.01 


0.65 




ngc2699 


4 


21.8 fl 


<0.1 












ngc2768 


4 


16.7 T 


10.71 ±0.02 


09 11 37.418 


60 02 14.84 


3.59 ± 0.01 


0.54 


14.5 ± 0.6 


ngc2778 





25.4 T 


<0.1 












ngc2974 


3 


25.9 T 


5.22 ± 0.02 


09 42 33.310 


-03 41 57.09 


4.19 ±0.02 


0.93 


10.4 ± 0.5 


ngc2986* 





22.3 T 


8.40 ± 0.03 


09 44 27.256 


-21 16 11.23 




160.14 




ngc3078 


4 


29.0 R 


124.95 ± 0.04 


09 58 24.630 


-26 55 36.09 


125.73 ± 0.04 


1.45 


279 ± 8 


ES0437-15 


3 


32.3 fl 


1.76 ±0.04 


10 36 58.100 


-28 10 34.70 


2.20 ± 0.05 


0.80 


3.2 ± 0.6 


ngc3156 


2 


14.0 T 


<0.1 












ngc3226 


3 


17.3 H 


7.29 ± 0.05 


10 23 27.005 


19 53 54.75 


2.61 ± 0.02 


0.97 




ngc3348 





38.5 fl 


1.66 ±0.02 


10 47 10.000 


72 50 22.71 


2.94 ± 0.04 


1.36 


7.8 ±0.5 


ngc3377 


1 


9.1 T 


<0.1 












ESO378-20 





35.6 R 


<0.1 












ngc3595 





30.4 R 


0.22 ± 0.01 


11 15 25.180 


47 26 50.60 


0.24 ± 0.01 


3.87 




ngc3610* 


3 


26.8 fl 


1.17 ±0.03 


11 18 20.700 


58 49 38.11 




230.78 




ngc4125* 


4 


20. 1 H 


1.23 ± 0.02 


12 08 04.180 


65 09 41.32 




86.29 




ngc4233 


4 


29.6 R 


2.52 ± 0.01 


12 17 07.679 


07 37 27.33 


2.64 ± 0.01 


1.02 


2.9 ± 0.5 


ngc4365 





15.7 R 


<0.1 












ngc4406 


4 


n.o v 


0.59 ± 0.02 


12 26 11.770 


12 56 46.40 


0.204 ± 0.07 


1.37 




ngc4476 


3 


24.7 T 


<0.1 












ngc4494 


4 


17.8 H 


0.27 ± 0.01 


12 31 24.030 


25 46 30.01 


0.10 ±0.01 


2.00 




ngc4552 


4 


n.o v 


93.40 ± 0.02 


12 35 39.805 


12 33 22.78 


32.30 ± 0.01 


0.35 


100 ±3 


ngc4697 


4 


15.5 T 


<0.1 












ngc4742 


3 


15.9 T 


<0.1 












ngc5198 





34.1^ 


0.83 ± 0.02 


13 30 11.390 


46 40 14.80 


1.15 ± 0.03 


1.16 


3.6 ± 0.4 


ngc5322 


4 


23.9 T 


13.60 ± 0.02 


13 49 15.269 


60 11 25.92 


9.33 ± 0.01 


1.08 


64 ±2 


ngc5557 





42.5 H 


<0.1 












ngc5576 





19.1 H 


<0.1 












ngc5812 


4 


24.6 fl 


<0.1 












ngc5813 


4 


24.6 fl 


2.95 ± 0.02 


15 01 11.234 


01 42 07.10 


2.14 ±0.01 


0.72 


12.3 ± 0.7 


ngc5845 


4 


18.1 T 


<0.1 












ngc5982 





39.3 fl 


<0.1 












ngc6278 





37.1 R 


1.06 ±0.01 


17 00 50.325 


23 00 39.73 


1.75 ± 0.02 


0.62 





Col.(l): name of the galaxy, sources with a star are far from the nuclear region of the corresponding galaxies and were treated as non detections; 
Col. (2): level of dust: = no dust, 1 = filamentary low, 2 = filamentary medium, 3 = filamentary high, 4 = dusty disk (Tran et al. 2001); Col. (3): 
distance in Mpc from (P) - Perrett et al. 1997, (T) - Tran et al. 2001, (R) - Rest et al. 2001, (V) - Virgo galaxies, assumed to be at distance of 
17 Mpc; Col. (4): flux at 3.6cm in mjy, or 4a upper limits for nondetections; Cols. (5) and (6): radio position (h, m, s) and (deg, arcmin, arcsec) 
from our maps (J2000); Col. (7): luminosity in 10 20 WHz -1 ; Col. (8): offset in arcseconds, between the 3.6 cm radio position and the position 
of the galaxy optical nucleus on HST images (Tran et al. 2001, Rest et al. 2001); Col.(9): peak flux from NVSS survey (Condon et al. 1998) 



non-AGN mechanism, e.g. nuclear starbursts. Since we are in- 
terested in the AGN/dust connection we wish to exclude this 
possibility. We argue here that the dominant source of radio 
emission in our detections is a non-thermal mechanism similar 
to that which operates in more powerful radio sources. 

There are several radio and infrared criteria that can be used 
to distinguish between emission from starburst and AGN galax- 
ies: (i) radio morphology, (ii) far-infrared to radio flux-density 
parameter u = log(56o^ m /<S'i.4GHz), (iii) infrared spectral 
index a IR = log^f^m/S^m)/ log(60/25) (Condon & 
Broderick 1988 and 1991; Condon et al. 1991; Condon, Huang 
& Thuan, 1991), and (iv) the steepness of the radio spectra. 



Radio morphology implies coherent radio jets and radio lobes 
that may lie well outside the optical galaxy. Starburst galax- 
ies usually have u > 1.6, and am > ±1.25. Steepness of 
the radio spectrum is also used as a criterion since optically 
thick AGN cores usually have fiat spectra, while the dominant 
emission from star-forming regions (supernova remnants, and 
cosmic rays diffusing from them) have steep spectra. Nearly 
all spirals and unclassifiable objects (e.g mergers) have steep 
spectra, while flat spectra and other AGN characteristics (radio 
morphology, u < 1.6, and ccir < ±1.25) are associated with 
ellipticals (Sadler, Jenkins, & Kotanyi 1989, Condon 1991). 
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All detected galaxies in our sample have low-luminosity 
unresolved sources in the innermost central regions. Although 
the sources are certainly nuclear in origin (suggesting AGN ac- 
tivity) any radio classification according to radio morphology 
is not possible (except in the clear case of a jet in NGC 5233). 
Half of the galaxies were picked based on their large scale 
dust and infrared properties from IRAS survey. This means 
that those galaxies are going to have larger am indices, which 
would mark them as starburst, although they still might have 
nuclear AGN which are the subject under discussion. Using 
the large scale IR emission to determine the nature of the nu- 
clear radio emission doesn't seem to be a very good discrim- 
inator between SBs and AGNs. However, most of our objects 
have measured nuclear Ha fluxes or upper limits (Tran, private 
communication), and standard calculations (Osterbrock 1989) 
indicate that the free-free fluxes from these regions would be 
below 3 /uJy, which is about two to three orders of magnitude 
smaller than our observed fluxes. Other evidence that we are 
dealing with non-thermal radiation comes from the flatness of 
the spectra in our sample. Eleven of the galaxies were detected 
before in the NVSS (Condon et al. 1998) and comparing the 
fluxes at our frequency (8.45 GHz) and the frequency of the 
NVSS (1.4 GHz) it is clear that most of the detected galaxies 
have flat spectra (Table 1). 

Previous studies (Phillips et al. 1986; Sadler, Jenkins & 
Kotanyi 1989) have shown that HII regions in early type galax- 
ies are not likely to contribute to the radio galaxy population 
above 10 19 W Hz -1 . Keeping in mind that all galaxies in our 
sample are Es and SOs, that emission is confined to nuclei of 
the host galaxies, and that the sources have flat spectra, we 
can assume that the dominant radio component in our case is 
synchrotron emission from an active nucleus producing low- 
luminosity counterpart of more distant, luminous AGNs. 

4.2. Radio Luminosity Function 

The size of our sample is too small and too limited in radio lu- 
minosity range to construct a complete local radio luminosity 
function (RLF) of early-type galaxies. In any case, the sam- 
ple was not constructed for that purpose. Still, we can make 
a useful estimate of the low luminosity end of the local RLF 
in order to see how it corresponds with previously found local 
RLFs and offer an estimate of the behavior of RLF at low lu- 
minosities. For this purpose we define the fractional luminosity 
function (Auriemma et al. 1977): 

Fi{L,z) = pi(L,z)/ipi(z), (1) 

where <pi(z) is the volume density of objects of type i at the 
redshift z, and pi(L,z) is the density of sources associated 
with optical objects of type i with the given radio luminos- 
ity L and at the given redshift z. The fraction of all ellipti- 
cal galaxies with luminosity at given frequency between L and 
L + dL at the redshift z, of optical magnitude M is given then 
by Fe,m{L, z)dL. In order to estimate the "bivariate" RLF 
defined in this way we can calculate the fractional detection 
fij = nij/Nij, where n(Lj, Mj) is the number of actually de- 
tected galaxies within the optical magnitude range Mj ± 0.5 
and radio luminosity interval logLj ±0.2, while N(L l , Mj) is 
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Fig. 1. The integral luminosity function derived from 17 de- 
tected sources out of a total of 36. The dots represent a crude 
LF calculated from the detections as the integral of a series of 
delta functions. A detection at L, contributes 1/Nd(i») where 
the denominator is the number surveyed galaxies detectable at 
Lj. The steep rise flattens off below of 10 20 W Hz -1 . The error 
bars are not plotted since the bins in the integral RLF are not 
independent. 

the number of galaxies in the sample which could have been 
detected if their optical magnitude and radio luminosity were 
in the given interval. In our case we did not bin in optical mag- 
nitude but only in radio power since the sample is limited. Our 
first estimate of the RLF is then given by Jl = hl /Nl and it 
is shown in Figure 1 in integral form, F(>L). As it can be seen 
from the Figure 1, the integral RLF of nearby ellipticals rises 
steeply with decreasing of the radio luminosity and only at the 
lowest intensities (10 19 W Hz -1 ) levels off at a point where 
<~ 50% of all E/S0 galaxies show activity. 

Previous RLFs (Auriemma et al. 1977; Sadler, Jenkins, & 
Kotanyi 1989; Condon 1991; Sadler et al. 2002) of nearby el- 
lipticals with AGN signature were made for galaxies with radio 
luminosities higher than 10 21 — 10 22 W Hz -1 . The more recent 
studies considered also starburst galaxies. While AGN were 
found in ellipticals, starbursts inhabited spirals. These differ- 
ent distributions had different RLFs and often starburst RLFs 
extended to the level of 10 20 W Hz -1 . With our low luminosity 
data, we are able to extend the existing RLFs of AGN down to 
10 19 W Hz -1 and can construct an exclusively AGN RLF. 

We compare our data with two studies (Sadler, Jenkins, & 
Kotanyi 1989 and Condon 1991) in Figure 2 (AGNs and star- 
bursts plotted). We have converted our differential data from 
F{L) to a spatial density <fi (number of sources per Mpc 3 per 
0.4 in log L) using the value for spatial density of early type 
galaxies, from Sadler, Jenkins, & Kotanyi (1989), which is 
jq-2.33 mag- 1 Mpc -3 . Gratifyingly our data agrees quite well 
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Fig. 2. Comparison of AGN and starburst (SB) local radio lu- 
minosity function. Filled symbols (circles: AGN, triangles: SB) 
are data from Sadler et al. (2002) and Condon (1991), while 
open circles are our data. The local density of AGN rises con- 
tinuously at low luminosities, reaching the value of SB, sug- 
gesting that AGNs are as common as SB in local universe. It 
is possible that at this low luminosity level both processes are 
present in galaxies, but in some galaxies one of the engines is 
stronger. 



Fig. 3. Plot of log radio luminosity in W Hz -1 versus abso- 
lute optical magnitude. Filled symbols are radio-detected nu- 
clei of galaxies, while open symbols indicate upper limits for 
the rest of the galaxies. Triangles are sources in galaxies with- 
out dust and circles are sources in dusty galaxies. Sources in 
dusty galaxies have a slight tendency for being more powerful 
than the sources in non-dusty galaxies. Luminosity error bars 
are smaller than the symbols 



with the previous data in the region of overlap. Together these 
data confirm the flattening of the RLF for AGNs below ~ 10 20 
W Hz -1 . It is also interesting that the space density of low lu- 
minosity AGN is very similar to starbursts galaxies of the same 
luminosities. The RLFs of the two distributions are basically 
overlapping in this luminosity range. 

5. Correlation of dust with radio emission 

5. 1. Crude Statistics 

The HST pictures of the galaxies in the original sample (Tran 
et al. 2001) confirm that dust is very common in ellipticals. 
There are two different morphologies in which dust appears in 
the galaxies from our sample: disky and filamentary. We have 
15 galaxies with disks and 9 with large amount of dust in fila- 
ments. Thirteen of the 24 dusty galaxies have a radio detection 
(54%), while 4 out of 12 non-dusty galaxies show a detection 
(33%). There is no significant difference in radio luminosity 
between the galaxies with disky and filamentary dust: 60% de- 
tections in galaxies with disks and 44% in galaxies with fila- 
ments. This finding is in general agreement with the findings 
by Tran et al. (2001). 

The relationship between optical absolute magnitude and 
radio luminosity for our weak radio sources is shown in Figure 
3. There is little difference in the distributions of the dusty and 



non-dusty galaxies, except perhaps that the three most powerful 
galaxies are all dusty. As expected, the more powerful radio 
sources are found in the brighter galaxies. 

Most nearby high luminosity radio sources are found in 
dusty early type galaxies (de Koff et al.(2000)), which suggests 
a link between dust and the existence of a radio source. Our 
wish now is to see if at the lower levels of radio luminosity 
dust also plays an important role. We divide the detections in 
two sets of sources: dusty and non-dusty according to the de- 
scriptions in Tran et al. (2001). As we see above, the dusty 
galaxies show a somewhat higher detection rate, but, given the 
steepness of the RLF, this could be influenced by slight differ- 
ences in the distances to the two samples, or slight differences 
in the achieved sensitivities. Therefore it is more meaningful to 
compare the RLFs of the two samples than the detection per- 
centages. 

5.2. Comparison of RLFs 

The integral RLFs for the two samples, computed by the same 
algorithm as that in Fig. 1 for the whole sample, are shown in 
Figure 4. In this representation also, the dusty galaxies seem 
more active, but the difference is relatively small (a factor of 
~ 1.6) and we wish to test the significance of this difference. 

We have tried two statistical tests: Kolmogorov-Smirnoff 
(K-S) and a test using maximum likelihood method (ML). The 
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Table 2. Estimated integral luminosty coefficients. Values of coefficeints of estemated integral luminosity function F obtained by 
maximum likelihood method. Errors are la estimates. 



data sample 

all data 
dust data 
non-dust data 



ct 



— 0.6 ± 0.3 
— 0.6 ± 0.4 
-0.5 ±0.3 



0.80 ±0.10 
0.74 ±0.09 
1.20 ±0.60 



20.41 ± 0.04 
20.40 ±0.10 
20.24 ± 0.04 



A 



0.21 ±0.04 
0.25 ±0.07 
0.14 ±0.02 




Fig. 4. The separated integral luminosity function. Open cir- 
cles present the RLF for sources from galaxies lacking dust. 
Filled circles present the RLF of sources from dusty galaxies. 
Statistical tests show that the two distributions are not distin- 
guishable, suggesting that dust is not important for the exis- 
tence of low-luminosity AGN in nearby early-type galaxies 



K-S test has the advantage of being parameter and form free, 
but the disadvantage of not being very conclusive for small 
samples. We have two data sets, one with 13 sources in dusty 
galaxies and one with 4 sources in non-dusty galaxies. We used 
routines from Numerical Recipes (Press et al. 1992). The prob- 
ability that these two observation sets could be obtained from 
the same RLF is 64%, hence the RLFs are statistically indis- 
tinguishable. However, the K-S test is sensitive to the effective 
number of data points, N e , which in our case of two distribu- 
tions is N e = NiN 2 /(Ni + N 2 ) = 3.1. Press et al. (1992) give 
N e > 4 as a limit for a decent accuracy. Thus, the above prob- 
ability is not very accurate, but it still implies that the two data 
sets (two luminosity functions) are not significantly different. 

Another approach that is more sensitive, but requires more 
a priori assumptions, is to fit a specific, parameterized, function 
to the RLF data using the maximum likelihood method, and 
compare the fitted functions. Since there are a limited number 
of degrees of freedom in the function, more powerful statistical 



Fig. 5. Estimated integral luminosity functions compared to 
the data. The drawn lines show the ML model fits (equations 
2 and 3) to the data, while the points are computed from in- 
dividual detections. Field circles present dusty galaxies while 
open circles present non-dusty galaxies. The thick line shows 
the ML model fit for dusty galaxies and the thin line for non- 
dusty galaxies. 

statements can be made. To estimate the integral luminosity 
function we used a set of two power-law functions allowing for 
a break in the RLF. Our choice is similar to some previously 
used functions (Auriemma et al. 1977): 

F = A ■ lO-P(x-Xc) for x > x c (2) 

= A)-(l + -(10" a(x " Xc) -l)) forx<x c (3) 

a 

where x = logi L, and L is radio luminosity in W Hz 1 . The 
normalization constant A is chosen so that at x = x c , F = A . 
Originally we assumed that at low luminosities the value of 
F(x) had to approach F = 1 as x — > — oo, thus providing an 
additional constraint on the model. These solutions provided 
poor fits to the data and were dropped. This implies, however, 
that there is another break (or continuous change of slope) in 
the RLFs below the limits of our survey. 
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The system of coefficients a, (5 (slopes of the curve), x c (the 
position of the break) and A (normalization) that maximize 
the probability in the method, provide also the best fit to the 
data. Table 2 contains the calculated values for a, (3, x c and 
A . 

The best-fit integral luminosity functions are compared to 
the observed values in Figure 5. The symbols are filled circles 
for dusty and open circles for non-dusty sources. The thick 
represents the model fit to dusty sources, while the thin line 
shows the fit to the non-dusty. The model curves fit the indi- 
vidual points reasonably well. The three radio brightest galax- 
ies lie somewhat above the best two power-law fit in the re- 
gion of the break, but the best fit value of the slope above the 
break, (3 ~ 0.75, agrees with the slope measured by Sadler 
et al. (2002) (Fig. 2) based on much more data in the higher 
luminosity ranges. 

The ML parameters a, (3, and x c are essentially identical 
for dusty and non-dusty galaxies, indicating that the forms of 
the RLF are similar. Not surprisingly the normalization Ao is 
higher for the dusty galaxies by a factor of about 1 .8, but this 
is only 1.6 times the uncertainty. 

Another way to globally judge the significance of the differ- 
ence between these RLFs it to ask if the "true" RLF were given 
by the dusty model, how unlikely is it that we would only detect 
four (or less) of the twelve non-dusty galaxies. If this probabil- 
ity is small, then the samples are significantly different. From 
Poisson statistics the probability of 4 or less non-dusty detec- 
tions given the dusty RLF (thick line on Fig. 5 or the second 
line values in the Table 2) is 27%, indicating a low statisti- 
cal significance. If, hypothetically, the true non-dusty RLF is a 
factor of 1.6 lower than the dusty RLF we can ask how many 
non-dusty galaxies must be surveyed in order to demonstrate 
the RLF difference at a reliability of, say, 5%. Repeating the 
Poisson analysis indicates that a sample about four times the 
current size is needed, or about 50 non-dusty galaxies. 

Perhaps the most important result of this investigation is 
that in any case, a sizable fraction of the non-dusty galaxies, 
~ 30%, are radio emitters, so that the presence of visible dust 
is not necessary for radio emission from an AGN. 

6. Conclusions 

We report 3.6 cm VLA observations of a sample of 36 near- 
by ellipticals selected on their optical/IR properties. We de- 
tected 17 unresolved (except the jet in NGC5322), compact, 
flat-spectrum radio cores associated with the central 1" of the 
nuclei, suggesting that all detected sources are low luminosity 
AGNs. The lowest detected luminosities are ~ 10 19 W Hz -1 . 

We determine the Radio Luminosity Function (RLF) from 
these galaxies down to a luminosity almost two orders of mag- 
nitude lower in luminosity than previously published studies. It 
shows the continuation in the rise of space density of sources 
with AGN signature, which was expected from other, unpub- 
lished, studies (Condon, private communication). At the lumi- 
nosities considered (i.e L ~ 10 19 - 10 22 WHz -1 ), the space 
densities of the AGNs and starburst galaxies approach each 
other, becoming hardly distinguishable. At the lower luminos- 



ity end of our sample ~ 50% of E/S0 galaxies have detectable 
radio-AGNs. 

Although the non-dusty galaxies show an indication of a 
lower probability of radio emission, the difference is not sta- 
tistically significant in a sample of this size. Dust detectable in 
HST images is certainly not necessary for nuclear radio emis- 
sion. This situation may be different for the more powerful ra- 
dio galaxies observed in earlier surveys. 

This takes us back to the question of fuel for the central 
engine of our low luminosity AGNs. If fuel is necessary for 
nuclear activity why do we find weak AGNs without visible 
dust? It should be noted that our non-dusty galaxies with ra- 
dio detections lay further away (D > 30 Mpc) and it might be 
possible that small amounts of dust were not detected. Also, ex- 
tremely diffuse gas and dust would not be visible (Goudfrooij 
& de Jong, 1995) but current theories of accretion require bars 
and disks and other distinct structures, so fueling from diffuse 
gas seems unlikely. Similarly, these galaxies could be fueled 
by hot, dust-free gas, but it seems unlikely that any mechanism 
in these low-luminosity sources would destroy dust more than 
in the high luminosity sources were dust is common. A more 
likely explanation is that the amount of dust and gas present 
near the nucleus is in some sense positively correlated to the 
AGN luminosity. The sources in our study are two to three or- 
ders of magnitude less luminous than the 3C sources in de Koff 
et al. (2000), where typical dust optical depths were of order 
unity. In the HST images, optical depths of less than 1 % would 
probably be missed. Alternatively, AGN fueling may be cyclic, 
and AGN radio emission is now fueled from material at a few 
Schwarzschild radii, after the material in a larger circumnuclear 
accretion disk has been temporarily consumed. 

The luminosity of an AGN is determined by the fueling 
rate and the mass-to-radiation conversion efficiency. The latter 
is influenced by the degree of advection which is in turn in- 
fluenced by the Eddington luminosity and the mass of the BH. 
As recent evidence suggests (Ho, 2002 and references therein) 
many of the characteristics of low luminosity AGN could be 
explained by an advection dominated accretion flow (Narayan 
& Yi 1995; Narayan, Mahadevan, & Quataert 1998). The ex- 
planation of the dust/radio emission/luminosity relations may 
perhaps be found when we know the BH masses of the galax- 
ies, or when we understand the characteristics of non-steady 
accretion flows. 
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